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Abstract 

In spite of its rather weird properties which include violation of the dominant-energy 
condition, the requirement of superluminal sound speed and increasing vacuum-energy 
density, phantom energy has recently attracted a lot of scientific and popular interests. 
In this letter it is shown that in the framework of a general k-essence model, vacuum- 
phantom energy leads to a cosmological scenario having negative sound speed and a 
big-rip singularity, where the field potential also blows up, which might occur at an 
almost arbitrarily near time in the future that can still be comfortably accommodated 
within current observational constraints. 
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The existence of phantom dark energy in the universe actually is a possibility not 
excluded by observations which has recently been widely discussed [1]. The physical 
properties of vacuum phantom energy are rather weird, as they include violation of the 
dominant-energy condition, P + p < 0,, naive superluminal sound speed and increasing 
vacuum-energy density. The latter property ultimately leads to the emergence of a 
singularity - usually referred to as big rip - in a finite time in future where both the 
scale factor and the vacuum-energy density blow up [2] . The existence of a singularity 
in finite time was already considered by Barrow, Galloway and Tipler in 1986 [3], even 
under the much weaker conditions p > and p + 3P > 0, by assuming that dP/dp 
is not a continuous function. This can actually be regarded as the first example of 
a big rip singularity. On the other hand, if we want the weak energy condition to 
be preserved one must regard the stuff of phantom energy to be made up of axions, 
at least when dealing with a quintessence field [4]. Indeed, if a quintessential scalar 
field with constant equation of state P^ = up^ is considered, then phantom energy 
can be introduced by allowing violation of dominant energy condition, P^ + p^ < 0, 
or what is equivalent, rotation of to imaginary values, — > i<£>, in the Lorentzian 
manifold [Notice that if the pressure and the energy density are respectively defined 
as P<f, = \4> 2 — V{4>) and p^ = |0 2 + V{4>), with V(4>) the potential energy, then it 
follows that Pcf, + p^ = (1 + Lo)p^, = 2 , and hence 2 < (i.e. classically an axionic 
component for vacuum phantom energy [4]) if we want the weak energy condition to 
be satisfied also for u; < —1, which will in turn automatically ensure violation of the 
dominant-energy condition]. No extra constraints are imposed to ensuring a causal 
propagation condition that dPj dp does not exceed unity. 

I will argue however that whereas P^ + p^ < or — > i$ will suffice to ensure 
uj < — 1 and phantom energy, or vice versa, in all quintessence cases studied so far, the 
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eventual emergence of a big rip in the future will take only place in scalar-field models 
with equations of state of the simplest form P^ = up^. In fact, in case of a Chaplygin 
gas with equation of state P^ = —A/p 1 ^ (where A and n are positive constants), which 
can in fact be regarded as just a particular case of a dust fluid with a given bulk 
viscosity in a k=0 Friedmann- Robertson- Walker universe [5] , the existence of phantom 
energy does not lead to emergence of a big rip [6]. 

In this short letter I will show that a cosmological model with a singular big rip 
at an arbitrary finite time in the future can be also obtained when the scalar field 
satisfies equivalent phantom-energy conditions in the case that it is equipped with 
non-canonical kinetic energy for models restricted by a Lagrangian of the form 

L = K(<P)q(x), (1) 

where x = ^A^A^. Such a definition, which of course includes the quintessence 
model as a limiting case, generally describes more general models claimed to solve 
the coincidence problem without fine tuning, which have been dubbed as k-essence 
[7]. Some of the current k-essence models featured suitable tracking behaviour during 
radiation domination with further attractors [7]. Introducing the usual variable y = 
1/y/x and re-expressing q(x) as q[x(y)] = g(y)/y, from the perfect-fluid analogy, we 
have for the pressure and energy density of a generic k-essence scalar field [8] 

PM - ^ (2) 
<>M = -K(M(y), (3) 

where the prime means derivative with respect to y. Now, the equation of state param- 
eter and the effective sound speed can be shown to be given by the i^(0)-independent 
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expressions 



cUy) 



_ P<t> _ g(y) 
P4> V9'{y) 
p± = g{y) - yg'jy) 
p'4, y 2 g"{y) 



(4) 



(5) 



In general, k-essence models are defined by taking K(4>) = 0~ 2 > [8]. Thus, for 
the weak energy condition to hold it follows from Eq. (3) that in these models the 
function g(y) must be decreasing. Moreover, in these models it is currently assumed 
that c 2 S( p > and hence Eq. (5) implies that g"(y) > 0, i.e. g(y) should be a decreasing 
convex function [8]. 

We set next the general form of the function g(y) when we consider a phantom- 
energy k-essence field; i.e. when we introduce the following two phantom-energy con- 
ditions: K(4>) < and 



which are just the conditions that would follow, both at once, whenever is made 
imaginary as in the quintessence models [4]. However, since the kinetic term is non- 
canonical in the k-essence scenario, the above two conditions should be defined by 
themselves, not as being derived from the general formalism of k-essence by Wick 
rotating the scalar field, for otherwise both the variable y and hence the function g(y) 
would turn out to be no longer real. In what follows we shall therefore introduce the 
above two phantom-energy conditions while keeping y and g(y) real. 

Now, the first of these conditions and Eq. (3) amount to g'(y) > in order for 
satisfying the weak energy condition > 0, and then from g'(y) > and the second 
phantom-energy condition, we deduce that g(y) > yg'{y). Whence using g'{y) > 0, it 
also follows that g(y) > 0. Therefore the function g(y) should be an increasing concave 
function, that is we must also set g"(y) < 0. We have then from Eq. (5) that the 



P<t>(y) + pM = 2K((f>)xdq(x)/dx < 0, 
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square of the speed of sound should necessarily be definite negative. However, even 
though an imaginary sound speed would at first sight mean catastrophic accelerated 
collapse of inhomogeneities, such a kind of instability could still be avoided at least at 
the subhorizon scale by taking into account the dependence of the sound speed on the 
wavelength characterizing the instabilities [9]. Finally, it is also a consequence from 
the above two phantom-energy conditions that oJ^iy) < —1. 

A simplest family of ^-functions satisfying the above requirements is (see Fig. 1 (I)) 



with B and (3 being given constants such that B > and < (3 < 1. Actually a 
more general function g(y) can be written as a polynomials g(y) = J2iBi% / S where 
the first term is given by Eq. (6) and all other extra terms are characterized with 
powers < $ < 1 as well and coefficients B > B\ > B 2 > .... It is moreover worth 
mentioning that the polynomial g(y) cases seem to be linked to the eight asymptotes 
discussed by Barrow [10] when applying Fowler theorems for first-order differential 
equations to obtain solutions of the Raychaudhuri equation which are continuous, finite 
and monotonic as t — > oo. Even though rigorously checking whether or not such a 
connection actually exist is outside the scope of the present work, it would appear 
interesting to investigate it. However, for the aims of this letter it will suffice taking 
only the first term of such a polynomials. Let us specialize then in the case of a spatially 
flat Friedmann- Robert son- Walker spacetime with line element 



in which a(t) is the scale factor. In the case of a universe dominated by a k-essence 
phantom vacuum energy, the Einstein field equations are then 



g(y) = By?, 



(6) 



ds 2 = -dt 2 + a(t) 2 dr 2 , 



(7) 



3H 2 =p4y), 2H + P(P (y) + P^y) = 0, 



(8) 
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Figure 1: (I) Generic shape of the function g(y) for a k-essence vacuum phantom-energy field. All 
units in the plot are arbitrary. (II) Evolution of the scale factor a(t) with cosmological time t for a 
function g(y) with the form given in Fig. 1 (I). The dashed line for a universe with an ever decreasing 
size corresponds to a choice of the sign of the integration constant such that t* < 0. If we choose that 
sign to be positive (solid line), then the constant t* becomes an arbitrary time in the future at which 
the big rip takes place. All units in the plot of Fi§. 1 (II) are also arbitrary. 



with H — a/a, the overhead dot meaning time derivative, ' = d/dt. Combining the two 
expressions in Eq. (8) and using the equation of state we can obtain for the function 
g(y) as given by Eq. (6) 



a* = ^. (9) 



a0C T. . x2/3/[3(l-/3)] ' (10) 



For our spatially flat case we have then the solutions 

1 

(f-t. 

where t* is an arbitrary constant. If we choose t* < 0, then the scale factor would 
ever decrease with time (see Fig 1 (II), dashed line). Obviously this solution family 
does not represent an accelerating universe and should therefore be discarded. Of quite 
greater interest is the choice > for which the universe (Fig. 1 (II), solid line) will 
first expand to reach a big-rip singularity at the arbitrary time t — t* in the future, to 
thereafter steadily collapse to zero at infinity; that is it matches the behaviour expected 
for current quintessence models with w < — 1. The potentially dramatic difference is 
that whereas in quintessence models the time at which the big rip will occur depends 
nearly inversely on the absolute value of the state equation parameter, in the present 
k-essence model the time t* is a rather arbitrary parameter. 

In the case that we take for the field potential the usual expression K{<p) = 4>~ 2 [8], 
the Euler-Lagrange equation for the current k-essence field can also be written as 



2/ 3 ^^0 - 3Hy 



d 2 g(y) 



dy 



dg(y) 



-g(y) 
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dg{y) 



+ ^ = 0. (11) 

<p 



dy 

Therefore, using Eqs. (6) and (10) in the case that K(<j)) = <p~ 2 one can integrate Eq. 
(11) to obtain for the phantom-energy k-essence field 



13-1 



(12) 
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with D and E being arbitrary integration constants and a an also arbitrary pre- 
factor for the scale factor in Eq. (10). We notice that the phantom field tends to 
vanish as t — > i*, and hence its potential, V = K(<f>) = <fi~ 2 , blows up at the big rip, 
such as it happens in quintessence models [4]. It follows as well that the energy density 
for the phantom field will increase initially as t — > t* and blow up at t = t*, as one 
should expect. 

The main result in this letter is that phantom vacuum-energy leads to a big rip 
singularity also for k-essence dark energy. Moreover, one can play with the arbitrary 
values of the pre-factor ao for the scale factor expression in Eq. (10) and those un- 
boundedly small positive values of t* which satisfy the observational constraint [11] 
1 > (3 > 0.7 (note that in the present model (3 = and the currently observed 

cosmic acceleration rate [11] to check that such a set of present observations [11,12] is 
compatible with unboundedly small positive values of t*. Even though unboundedly 
larger values of t* are also allowed this way, k-essence phantom energy certainly may 
well allow a very near occurrence of the big rip in the future. Therefore, one could 
say that, in cosmo logical-time terms, a far and a near occurrence of the big rip are 
similarly probable and that, as a consequence from this, in the framework of cosmic 
k-essence, cosmic doomsday might be awaiting us around the corner. 
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